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Executive  Summary 


Impedance  matching  is  a  canonical  problem  in  electrical  engineering.  The 
problem  is  to  maximize  power  flowing  from  a  generator  to  a  load  with  an 
adroit  design  of  a  matching  network.  Typically,  matching  problems  use  a 
single-port  load.  In  contrast,  this  report  develops  matching  for  a  load  with 
multiple  ports.  The  specific  multiport  load  is  a  three-wire  antenna.  Each  wire 
has  a  feed  port  so  that  the  load  is  a  3-port.  The  matching  problem  is  to  And 
a  matching  network  that  maximizes  power  flow  from  a  single  feed  to  all  three 
ports  across  a  wide  frequency  band  (100-700  MHz).  Unexpectedly,  several 
theoretical  matching  multiports  were  found  that  gave  excellent  wideband 
performance.  Consequently,  this  serendipitous  result  opens  several  directions 
of  opportunity. 

One  direction  aims  at  enhancing  the  antenna’s  wideband  performance. 
Understanding  the  physical  phenomenon  of  coupling  between  the  wires,  the 
physical  geometry  of  the  wires,  and  a  Pareto  theory  for  maximal  bandwidth 
and  minimal  size  are  all  significant  research  topics.  Another  direction  aims 
at  the  matching  networks  because  of  the  gaps  in  theory  and  practical  design. 
Specific  and  detailed  formulation  of  these  research  opportunities  is  made 
explicit  in  this  report. 
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1  Why  Multiport  Matching? 

Impedance  matching  is  a  ubiquitous  problem  in  electrical  engineering.  The 
standard  matching  problem  searches  for  a  2-port  that  maximizes  the  power 
transfer  from  a  generator  to  a  load.  Increasing  bandwidth  in  communication 
systems,  the  integration  of  many  functions  on  a  chip,  and  multifunction 
antennas  all  lead  to  the  problem  of  multiport  matching.  In  general,  when 
a  generator  and  a  load  are  multiports,  the  connecting  matching  network 
must  also  be  a  multiport. 

Multiport  matching  assumes  the  ports  of  the  load  are  coupled.  Other¬ 
wise,  uncoupled  ports  reduces  matching  to  a  collection  of  standard  2-port 
matching  problems.  A  multiport  antenna  provides  a  splendid  case  study. 
When  antenna  arrays  are  spaced  closer  than  half  a  wavelength  (A/2),  mutual 
coupling  effects  emerge  [38] : 

When  using  multiple  antenna  elements  for  diversity  implemen¬ 
tation  on  small  personal  communications  devices,  the  resulting 
closely  spaced  antenna  elements  exhibit  the  well-known  mutual 
coupling  . . .  [that]  alters  both  their  terminal  impedance  and  ra¬ 
diation  pattern  characteristics.  These  changes  obviously  impact 
the  diversity  performance  of  the  multiantenna  system. 

A  powerful  design  approach  is  the  emerging  port  decoupling  [9] : 

[port  decoupling]  involves  a  modal  feed  network  which  makes  use 
of  the  orthogonality  of  the  eigenmodes  of  the  array  to  achieve 
decoupling.  The  input  ports  to  the  feed  network  and  array  com¬ 
bination  can  then  be  matched  independently. 

Figure  1  illustrates  this  decoupling  and  single-channel  matching  approach. 
Figure  2  shows  a  6-port  matching  circuit  where  the  decoupling  and  matching 
have  been  merged  into  a  nicely  interconnected  multiport.  Figure  3  shows 
a  matching  problem  where  a  3-port  antenna  is  powered  by  a  single  feed. 
The  4-port  matching  networks  will  exploit  the  mutual  coupling  to  attain  a 
wideband  antenna. 
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Figure  1:  Matching  and  decoupling  schematic  [37]  (c)  IEEE  2008. 


Figure  2:  6-port  decoupling  and  matching  [41]  (c)  IEEE  2006. 
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Figure  3:  4-port  matching  with  a  single  feed  at  Port  1. 
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Section  3  describes  the  three-wire  antenna,  the  three  feed  points,  and  the 
computation  of  its  3  x  3  admittance  matrix  over  100-700  MHz.  This  admit¬ 
tance  matrix  shows  the  resonances  of  the  wires  and  the  coupling  between  the 
wires. 

Section  4  starts  the  matching  process  by  directly  connecting  the  three  feed 
points  to  a  common  terminal.  Thus,  the  three-wire  antenna  is  turned  into  a 
standard  single-port  antenna  amenable  to  2-port  matching.  Unfortunately, 
this  direct  connection  has  a  terrible  Voltage  Standing  Wave  Ratio  (VSWR) 
over  100-700  MHz.  The  smallest  VSWR  attainable  by  any  matching  2- 
port — regardless  of  the  number  of  components  and  topology — exceeds  22. 
Only  17%  of  the  power  is  delivered  to  the  antenna — the  remaining  83%  is 
reflected  back  to  the  feed.  Multiport  matching  improves  this  power  delivery. 

Section  5  matches  of  the  three-wire  antenna  over  the  class  of  4-ports 
shown  in  Figure  3.  The  matching  results  are  excellent:  the  VSWR  is  less 
than  2.7  over  140-700  MHz.  This  multiport  matching  network  delivers  79% 
of  power  to  the  antenna;  only  21%  is  reflected  back  to  the  feed.  This  im¬ 
provement  of  this  wide-band  VSWR  motivated  this  exploration  of  multiport 
matching. 

Section  6  straddles  the  preceding  matching  topologies  by  matching  with 
a  “channelized”  4-port — each  port  of  the  antenna  is  connected  to  its  own 
matching  2-port.  The  preliminary  results  show  that  this  “semi-direct”  con¬ 
nection  is  not  as  viable,  and  the  cross-coupling  in  the  full  multiport  is  neces¬ 
sary  for  a  wideband  match.  This  observation  again  raises  questions  regarding 
a  physical  understanding  of  multiport  matching.  Section  7  organizes  these 
results  into  a  collection  of  research  topics: 

•  Developing  H°°  bounds  to  compute  absolute  matching  performance 

•  Extracting  a  circuit  from  the  general  multiports 

•  Considering  ladder-like  multiport  topologies 

•  Understanding  the  role  of  coupling  and  matching 

•  Other  multiport  matching  applications 

Appendix  A  offers  an  in-depth  description  of  the  resonances  of  the  three- 
wire  using  characteristic  modes.  Appendix  B  reviews  the  necessary  back¬ 
ground  on  the  lumped,  lossless  multiports  used  for  matching. 
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2  Notation 


Table  1  lists  common  notation.  The  last  three  entries  are  the  scattering, 
impedance,  and  admittance  matrices.  The  implicit  assumptions  when  using 
these  matrices  is  the  existence  and  linear  relationship  of  the  currents  and 
voltages  [28].  For  example,  the  admittance  matrix  Y  of  a  three-port  links 
the  current  and  voltage  vectors  as 
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1/21  1/22  1/23 
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CO 

bO 

CO 

CO 

Vs 

Table  1:  Notation. 


symbol 

Description 

M 

Real  numbers 

C 

complex  plane 

c+ 

open  right  half  complex  plane 

j 

positive  square  root  of  -1 

P 

complex  frequency:  p  =  a  +  ju 

a 

neper  frequency  (nepers  per  second) 

u 

radial  frequency  (radians  per  second) 

In 

N  X  N  identity  matrix 

S 

scattering  matrix  (unitless) 

Y 

admittance  matrix  (Siemens) 

Z 

impedance  matrix  (ohms) 

Figure  4  illustrates  the  cascade  of  a  2-port  with  scattering  matrix 

^  _  -Sii  Si2 

-S21  <§22 

with  Port  2  terminated  in  a  load  with  reflectance  sl-  The  reflectance  si 
looking  into  Port  1  is 

Sl  =  d{S,  Sl)  :=  Sii  -F  Si2Sl(1  -  S22Sl)~^S21. 
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When  the  dependence  on  the  complex  frequency  p  must  be  made  explicit, 
the  notation  takes  the  form 


Slip)  =  sl;p) 

:=  Sub)  +  Si2b)sLb)(l  -  S22ip)sLip))~^S2iip). 


h 


Figure  4:  Cascade  of  a  2-port  and  load. 


Let  a  multiport  have  a  scattering  matrix  with  a  block  structure 


“S'!!  512 
^21  S22 


Let  5ii  be  an  M  X  M  scattering  matrix  that  corresponds  to  Ports  1,  2  , 

.  Let  S22  he  duN  xN  scattering  matrix  that  corresponds  to  the  remaining 
ports.  If  these  remaining  N  ports  are  terminated  in  an  iV-port  load  with 
scattering  matrix  5^,  the  scattering  matrix  looking  into  Ports  1,  2,  . . . ,  M  is 

:SiS,  Sl)  :=  5n  +  5i25l(/v  -  5225l)-'52i, 

where  In  denotes  the  N  x  N  identity  matrix. 
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3  The  Three-Wire  Antenna 


An  arbitrary  three-wire  antenna  was  chosen  to  investigate  the  value  of  mul¬ 
tiport  matching.  Figure  5  shows  that  the  antenna  consists  of  three  closely 
spaced  center-fed  vertical  dipoles  operating  over  100-700  MHz. 


Figure  5:  Center-fed  dipoles. 

The  dipole  lengths  and  resonant  frequencies  are  listed  in  Table  2.  Com¬ 
paring  the  wavelengths  to  the  antenna  distances  justihes  “closely  spaced”  as 
a  descriptor. 


Table  2:  Dipole  specihcations  (meters). 


Wire 

X 

y 

Length 

1 

0.00 

0.00 

0.50 

2 

0.05 

0.00 

0.25 

3 

0.00 

0.05 

1.00 

The  admittance  matrix  of  this  three-wire  antenna  provides  insight  into  its 
multiport  response.  Specihcally,  the  admittance  matrix  shows  the  resonant 
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and  mutual  coupling  frequencies.  It  suffices  to  define^  these  frequencies  from 
the  local  maxima  of  admittance’s  magnitude.  The  admittance  matrix  relates 
the  feed-point  voltages  to  the  feed-point  currents  measured  at  Ports  1,  2,  and 
3  of  Figure  5: 


ii 

yii  yi2  yis 

Vi 

i2 

= 

t/21  t/22  t/23 

V2 

is 

ysi  ys2  yss 

Vs 

The  feed  point  admittance  were  calculated  using  Version  14  of  the  Expert 
MININEC  Broadcast  Professional  [32].  This  software  is  an  advanced  engi¬ 
neering  tool  for  the  design  and  analysis  of  wire  antennas.  The  wire  antennas 
are  assumed  to  be  a  collection  of  arbitrary  thin,  straight  wires  in  free  space  or 
over  a  ground  plane.  A  method  of  moments  algorithm  solves  for  the  current 
distribution  on  the  wires.  The  solution  for  current  distribution  is  based  on 
the  numerical  solution  of  an  integral  equation  representation  of  the  electric 
helds  [33].  The  process  of  solution  begins  with  thin-wire  assumptions  [36]: 

W-1  Antennas  must  be  modeled  as  a  group  of  one  or  more  straight  wires. 

W-2  The  wire  radius  is  very  small  with  respect  to  the  wire  length. 

W-3  Parallel  wires  should  be  several  wire  diameters  or  more  apart. 

W-4  The  radius  is  very  small  with  respect  to  segment  lengths. 

Table  3  tabulates  the  dipole’s  lengths,  radii,  and  segment  lengths.  Compared 
to  the  wavelength  A  =  0.43  meters  at  700  MHz,  this  antenna  model  is  within 
MININEC  assumptions  W-1,  W-2,  W-3,  and  W-4. 


Table  3:  MININEC  wire  models  (millimeters). 


Wire 

Length 

Radius 

Segment 

Separation 

1 

500 

1 

25 

50 

2 

250 

1 

25 

50 

3 

1,000 

1 

25 

50 

^Subsequent  reports  will  include  a  multiport  generalization  of  resonance  and  anti¬ 
resonance  defined  on  the  Smith  chart. 
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The  electric  field  is  determined  by  the  vector  magnetic  potential  and  the 
scalar  electric  potential.  The  two  potentials  are  calculated  from  potential 
integrals  that  are  solutions  of  the  Helmholtz  vector  and  scalar  wave  equa¬ 
tions.  The  integrands  of  the  potential  integrals  are  the  wire  current  and 
wire  charge  distributions.  Because  wires  are  thin,  the  currents  are  axially 
directed  (i.e.,  no  circumferential  currents  on  the  wires).  The  current  and 
charge  are  linked  by  the  equation  of  continuity.  Expert  MININEC  Broadcast 
Professional  assumes  the  following  boundary  condition: 

TE-0  The  tangential  electric  field  at  the  surface  of  a  perfect  conductor  is 
zero. 

Because  the  wires  are  assumed  to  be  thin,  this  assumption  forces  the  total 
axial  electric  field  on  the  wire  to  zero.  The  three  sources  of  the  tangential 
electric  field  on  the  wire  are: 

TE-1  currents  and  charges  on  the  wires  and  on  nearby  wires, 

TE-2  incoming  waves  from  distance  or  nearby  radiators, 

TE-3  local  sources  of  electric  field  on  the  wire. 

The  local  sources  are  voltage  sources  or  current  sources  that  connect  to  the 
wires.  By  summing  the  tangential  electric  field  components  at  each  segment 
on  the  wire  antenna  and  enforcing  the  zero  total  value,  an  integral  represen¬ 
tation  for  the  currents  and  charges  is  obtained. 

The  electric  field  integral  equation  is  solved  in  Expert  MININEC  Broad¬ 
cast  Professional  by  the  method  of  moments.  Basis  functions  are  chosen  to 
represent  the  unknown  currents  (i.e.,  triangular  basis  functions).  Testing 
functions  are  chosen  to  enforce  the  integral  equation  on  the  surface  of  the 
wires.  These  basis  and  testing  functions  yield  a  matrix  approximation  of 
the  integral  equations.  If  this  matrix  is  inverted  and  multiplied  by  the  local 
sources  of  electric  field,  the  complex  magnitudes  of  the  current  basis  func¬ 
tions  are  derived.  All  antenna  performance  parameters  can  be  determined 
from  the  derived  current  distribution,  including  the  feed  point  admittance. 

Figure  6  plots  the  admittance  matrix  as  a  function  of  frequency.  Each 
(m,  n)  subplot  plots  the  magnitude  of  its  corresponding  normalized  admit¬ 
tance  \ym,n\-  The  antenna  resonances  are  prominent  features  of  the  port  ad¬ 
mittances  1/14,  1/2,2,  and  1/3,3  along  the  main  diagonal:  The  0.5-meter  dipole 
has  |//i,i|  showing  a  peak  near  300  MHz.  The  0.25-meter  dipole  has  |//2,2| 


peaking  near  600  MHz.  The  l-meter  dipole  has  |?/3^3|  peaking  near  150  MHz 
and  450  MHz. 


A, 

200  400  600 


MHz 


200  400  600 

MHz 


Figure  6:  Admittance  matrix — normalized  magnitudes. 


The  plot  also  shows  both  the  coupling  and  reciprocity  (Y  =  Y'^)  between 
the  dipoles.  For  example,  the  voltage  vi  across  Port  1  to  0.5-meter  dipole 
also  produces  currents  from  Port  2  (the  0.25-meter  dipole)  at  600  MHz  and 
Port  3  (the  1-meter  dipole)  at  300  MHz.  Likewise,  voltage  V2  across  Port  2  to 
the  0.25-meter  dipole  produces  currents  from  Port  1  (the  0.5-meter  dipole) 
and  Port  3  (the  1-meter  dipole)  at  600  MHz.  Finally,  the  voltage  v^,  across 
Port  3  feeds  to  the  1-meter  dipole  currents  from  Port  1  at  300  MHz  and  in 
Port  2  at  600  MHz.  Consequently,  the  admittance  plot  does  reveal  some  of 
the  resonances  and  mutual  couplings. 

Figure  7  plots  the  scattering  matrix  the  three-wire  antenna  as  a  function 
of  frequency  /.  For  /  <  200  MHz,  the  reflectances  are  plotted  in  black. 
For  /  >  200  MHz,  the  reflectances  are  plotted  in  blue.  The  Smith  chart 
reveals  the  resonances  of  the  input  ports:  A  resonant  frequency  is  dehned 
where  s{jujx)  crosses  the  real  axis  in  the  upward  direction.  An  anti-resonant 
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frequency  uj^  is  defined  where  s{ju^)  crosses  the  real  axis  in  the  downward 
direction  [43].  For  example,  S3^3{ju)  explains  that  two  large  spikes  at  150 
and  450  MHz  in  the  admittance  plot  are  resonant  frequencies  while  the  spike 
at  300  MHz  is  an  anti-resonant  frequency. 


Port  1 


e 

o 

Oh 


Port  2 


Port  3 


Figure  7:  Scattering  Matrix. 

Finally,  because  this  scattering  matrix  is  the  “object  under  discussion,” 
the  limitations  of  the  antenna  model  are  made  explicit.  First,  the  limitations 
of  MININEC  are  recognized  [36]  and  a  convergence  study  was  not  performed 
[33].  Second,  the  dipole’s  lengths  were  chosen  for  convenience.  A  more 
judicious  choice  of  dipole  lengths  will  produce  more  resonances  in  100-700 
MHz.  Third,  the  sensitivity  of  this  antenna  to  physical  perturbations  has 
not  been  studied.  Nevertheless,  the  wavelength  and  the  physical  dimensions 
argue  that  this  three-wire  antenna  is  “good  enough”  as  a  multiport  matching 
exemplar. 
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4  The  Direct  Connection 


This  section  starts  the  matching  process  using  the  direct  connection  as  a 
baseline  case.  Figure  8  is  a  schematic  of  the  direct  connection.  The  green 
box  is  the  matching  4-port  that  simply  wires  all  three  ports  of  the  antenna 
to  the  single  feed. 


Figure  8:  4-port  matching;  direct  connection. 


This  direct  connection  converts  the  3-port  antenna  into  a  1-port  antenna 
amenable  to  standard  matching  as  shown  in  Figure  9.  Consequently,  the 
performance  obtained  by  matching  this  direct  connection  benchmarks  the 
performance  attained  by  the  multiport  matching  networks. 


Figure  9:  2-port  matching  of  the  direct  connection. 


Figure  10  plots  the  transducer  power  gain  of  the  direct  connection: 


—  1 


where  Si{juj)  is  the  reflectance  looking  into  the  feed  to  the  antenna.  In 
addition  to  the  expected  resonances  near  150,  300,  and  600  MHz,  a  resonance 
exists  near  450  MHz.  Appendix  A  shows  that  this  resonance  is  on  the  long 
wire.  Figure  11  plots  the  reflectance  of  the  direct  connection  on  the  Smith 
chart.  The  number  and  location  of  these  loops  near  the  boundary  of  the 
Smith  chart  indicate  that  matching  is  difficult. 
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Direct  Connection 


Figure  10:  Transducer  power  gain  of  the  direct  connection. 


Three  Wires:  Direct  Connection 
HF  Wliip 

100  </<  690  MHz 


Figure  11:  Smith  chart  of  the  direct  connection. 
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Figure  12  shows  the  difficulty  of  matching  the  direct  connection.  The 
plot  shows  the  matching  performance  of  various  classes  of  2-ports  connecting 
the  feed  point  to  a  generator  as  shown  in  Figure  9.  The  vertical  axis  is 
the  VSWR.  The  horizonal  axis  is  the  degree  d  of  the  matching  2-port  and 
is  equivalent  to  the  number  of  lumped  elements  in  the  2-port.  The  curve 
labeled  “low-pass  ladders”  plots  the  VSWR  as  a  function  of  the  number  of 
elements  and  attains  a  VSWR  in  excess  of  25.  The  solid  patch  marks  the  H°° 
bound — the  smallest  VSWR  attainable  by  any  lossless  2-port.  This  VSWR 
exceeds  22  so  that  this  direct  connection  is  useless.  Between  these  two  plots 
lies  the  “state-space”  curves.  This  curve  plots  the  matching  performance 
available  from  all  lumped,  lossless  2-ports  of  degree  d.  This  matching  class 
is  reviewed  in  Appendix  B  and  is  the  basis  for  the  multiport  matching  in  the 
next  section. 


Direct  Connection  Matching 


degree  t/:  140</<690MHz 


Figure  12:  Matching  the  direct  connection. 
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5  Multiport  Matching 

Figure  13  is  a  schematic  showing  the  voltages,  currents,  and  incident  and 
reflected  waves  of  a  3-port. 


Figure  13:  Scattering  waves  for  a  3-port. 


The  incident  wave  [4,  Eq.  4.25a],  [5,  page  234] 


and  the  reflected  wave  [4,  Eq.  4.25b],  [5,  page  234] 


?V2- 


are  defined  with  respect  to  the  normalizing^  matrix 

ro,i  0  0 


d/2: 


Rq  — 


0 

0 


ro,2  0 

0  ro,3 


^Two  accessible  books  on  the  scattering  formalism  are  Baher  [4]  and  Balabanian  & 

Bickhart  [5].  Baher  omits  the  factor  of  1/2  but  carries  this  rescaling  into  the  power  defi- 

1/2 

nitions.  Most  other  books  use  the  power-wave  normalization  [14]:  a  =  '  {v-|-  Zoi}/2, 

where  the  normalizing  matrix  Zq  =  Rq-\-  jXg  is  diagonal  with  diagonal  resistance  Rq  >  0 
and  reactance  Xq. 
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This  document  always  normalizes  with  respect  to  tq  =  50  ohms: 

■  1  0  0  ■ 

Ro=ro  0  10  =  ro/3. 

0  0  1 

The  scattering  matrix  Sl{p)  maps  the  incident  signal  a(p)  to  the  reflected 
signal  b(p): 

hi  Sii  S12  <Si3  Oi 

b  =  62  =  >521  >522  >523  ^2  =  S^S.. 

_  hs  _  _  >531  S32  >533  _  _  0-3  _ 

The  complex  power^  delivered  to  the  3-port  is  [5,  page  241]: 

W{p)  :=  v(p)'^i(p). 

The  average  power  delivered  to  the  3-port  is  [18,  page  19],  [25]: 

Pl  ■=  ISIH'-]  =  -  b"b}  =  V{/  -  Sf  St}*-  (1) 

The  matrix  I  —  S^Sl  is  called  the  dissipation  matrix  [10]  or  the  radiation  ma¬ 
trix  [37].  The  latter  name  is  misleading  because  it  implies  that  the  power  Pl 
that  is  delivered  to  the  antenna  is  actually  radiated.  The  eigen-decomposition 
of  this  matrix,  which  is  equivalent  to  the  eigen-decomposition  of  the  scat¬ 
tering  matrix  Sl,  is  basic  to  matching  [37].  If  the  3-port  consumes  power 
(Pl  >  0)  for  all  its  voltage  and  current  pairs,  the  iV-port  is  said  to  be  passive. 
Because  v(p)  is  the  Fourier  transform  of  the  voltage  vector  V(t): 

/OO 

e-^'2-piv(t)df, 

■OO 

v(p)  has  units  V/Hz.  Likewise,  i(p)  has  units  A/Hz,  so  that  hF(p)  units  of 
W/Hz^.  In  comparison,  both  the  incident  wave  a  and  the  reflected  wave  b 
have  units  of  v^W /Hz  so  that  it  makes  sense  to  call 

Pa  =  a^a 

^Baher  uses  [4,  Eq.  2.17]:  W{p)  =  i{p)^v{p).  Gonzalez  [18]  power  definitions  require 
a  “1/2”. 
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the  incident  power  and 


n  =  b^b 

the  reflected  power,  respectively  [5,  page  227].  Optimal  feed  schemes  are 
those  that  minimize  the  reflected  power. 

min{Pb(a)  :  Pa  =  1}. 

Equivalently,  these  optimal  feed  schemes  are  also  solutions  to  the  matrix 
norm  problem: 

||/  -  SlSlW  =  max{a^(/  -  S'f  S'L)a  :  ||a||  =  1}. 

The  physical  interpretation  of  this  feed  scheme  is  that  it  maximizes  the  power 
delivered  to  the  antenna  subject  to  the  constraint  of  unit  incident  power  [37]. 
Figure  14  plots  the  norm  of  the  dissipation  matrix  as  a  function  of  frequency. 


ThreeWires  1  OOtoTOOMHz.:  Generalized  Transducer  Power  Gain 


Figure  14:  Unmatched  Gain. 
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The  matching  network  should  map  the  generator’s  power  to  an  output 
power  that  will  be  more  useful  at  the  load.  The  modihcation  of  power  is 
generically  called  “gain.”  The  matching  problem  requires  gain  computations, 
and  we  need  the  maximum  power  and  mismatch  dehnitions.  The  maximum 
power  available  from  a  generator  is  dehned  as  the  average  power  delivered 
by  the  generator  to  a  conjugately  matched  load  [18,  Eq.  2.6.7]: 

Pg, ma.x  ■=  ^^151=^  =  ^  2^  (1  ~ 

The  transducer  power  gain  is  [18,  page  213]: 

Pl  power  delivered  to  the  load 

(jx'jp  !—  -  —  - . 

Pg, raa.x  maximum  power  available  from  the  generator 

Lemma  1  Figure  15  is  a  schematic  of  a  f-port  that  matches  a  generator 
looking  into  Port  1  to  a  3-port  load  terminating  Ports  2,  3,  and  f. 


Figure  15:  Matching  4  port. 

The  reflectance  looking  into  Port  1  is 

Si  =  Sii  +  S'i25'l(J3  —  S'22*S'i)  ^S'21. 

The  transducer  power  gain  is 

Sg{h  -  S22SL)-^{h  -  S^SLflh  -  S22Sl)-^S21.  (2) 

|1  -  SgSiP 
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Proof:  The  first  task  is  to  link  the  power  Pl  delivered  to  the  3-port  load 
(Equation  1)  to  the  input  Oi.  The  4-port  scattering  equations  and  the  load 
reflectance  are: 


■  bi  ■ 

-Sll 

^12  ■ 

di 

b2 

^21 

822 

a2 

Sl^2  —  ^2- 


Eliminate  a2  from  the  lower  equation: 


b2  =  {h  —  82281)  ^82iai. 


The  power  Pl  delivered  to  the  3-port  load  is 

Pl  =  \^l8^l{h  -  8228l)-^{I3  -  ^l)(/3  -  ^22^L)“'^2iai. 

A  corresponding  symmetry  of  the  power  relations  exists  at  the  input  port. 
The  reflectances  at  Port  1  and  the  generator  are 


bi  =  siai]  ai  =  SGbi+CG. 


Eliminate  61  from  the  generator’s  equation: 


(1  -  SGSi)ai  =  cg- 

Substitute  into  the  power  Pl  to  the  load: 

Pl  =  -  8228L)-^{h  -  8^8L){h  -  8228l)-^82i. 

2  |1  -  SGSik 

Dividing  by  the  maximum  gain  available  from  the  generator  gives  the  result. 

/// 


Lemma  1  is  a  simple  algebraic  statement  that  holds  regardless  of  the 
passivity  of  the  matching  4-port  or  the  3-port  load.  All  that  is  required 
is  that  the  inverses  exist.  A  substantial  simplihcation  is  possible  when  the 
matching  4-port  is  lossless. 


18 


Corollary  1  Assume  the  schematic  of  Figure  15.  Assume  the  f-port  is  loss¬ 
less: 

S^S  =  h. 

Assume  the  load  Sl{p)  is  strictly  passive: 

SESl  <  h. 

Assume  the  reflectance  of  the  generator  is  zero: 

sg  =  0. 

The  transducer  power  gain  is 

Gt  =  1  —  |<sip.  (3) 

The  smallest  transducer  power  gain  over  a  frequency  band  hi  is  denoted 

||Gr(sG,  S',  S'l)||_oo,o  :=  min{Gr(sG,  S',  S'l;  jo;)  :  ct;  e  fi}, 

where  the  dependence  on  the  generator’s  reflectance  sg{p),  the  matching  4- 
port  S'(p),  and  the  load  S'l(p)  is  made  explicit.  The  generator’s  reflectance 
sg{p)  and  the  load  Sl{p)  are  fixed.  The  collection  of  matching  4-ports  speci- 
hed  by  the  circuit  designer  is  called  the  designable  part  and  denoted  ©.  The 
Single-Input  Multiple-Output  (SIMO)  matching  problem  is  to  find  lossless 
multiports  from  ©  that  maximize  the  transducer  power  gain  Gt  over  a  fre¬ 
quency  band  O.  Generalizing  this  SIMO  matching  problem  for  an  iV-port 
load  is  straightforward: 


SIMO  Matching:  Given  the  generator’s  reflectance  sg, 
an  7V-port  load  S'l(p)  and  a  class  ©  of  lossless  (iV  -|-  1)- 
ports,  maximize  the  transducer  power  gain; 

snp{\\GT{sG,S,SL)\\n,-o.:See}. 


A  general  matching  class  consists  of  all  the  lumped,  lossless  iV-ports  of 
degree  d.  This  class  is  denoted  by  d)  and  described  in  detail  in  Ap¬ 

pendix  B.  Physically,  G+(iV,  d)  represents  a// the  iV-ports  that  can  be  synthe¬ 
sized  with  at  most  d  lumped  elements — inductors  and  capacitors — connected 
only  by  wires,  transformers,  and  gyrators. 
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Mathematically,  U~^{N,  d)  consists  of  all  iV  x  iV  scattering  matrices  S{p) 
that  are  rational  fnnctions  analytic  and  bonnded  on  the  open  right  half  of 
the  complex  plane,  real  on  the  real  axis, 

S{p)  = 

nnitary  on  the  imaginary  axis 

S{juj)^  S{juj)  =  In, 
and  have  Smith-McMillan  degree 

degsM[-5(p)]  <  d. 

Figures  16,  17  18,  19,  20,  and  21  match  the  3-port  antenna  from  d). 
The  optimization  problem  is  the  SIMO  Matching: 

sup{||G't(sg,  ^L)|b,-oo  :  ^  e  f/+(4,  d)}. 

Figure  16  plots  the  4-port  matching  over  the  full  frequency  band  100-700 
MHz.  The  blue  curve  plots  the  unmatched  gain  from  Figure  14.  The  green 
curve  plots  the  gain  obtained  by  matching  over  77+  (4,  d)  when  d  =  A.  The 
transducer  power  gain  is  improved  from  -14  dB  to  -3  dB  but  the  poor  response 
of  the  antenna  near  100  MHz  limits  the  performance.  Figure  17  verihes  this 
statement  by  matching  from  the  d  =  A  multiports  but  limiting  the  frequency 
band  to  140-700  MHz.  Excellent  matching  performance  is  obtained  with 
the  nearly  flat  gain  function.  This  “flatness”  of  the  gain  at  optimum  is 
characteristic  of  theory. 

Conjecture:  SIMO  matching  solutions  are  characterized  by  flat 
gain  as  d  — cx). 

Figure  18  shows  that  matching  with  d  =  3  still  has  a  nearly  flat  response 
to  support  the  conjecture.  Figures  19,  20,  and  21  match  with  d  =  2,  1,  0 
with  only  a  0.6-dB  decrease  in  performance.  However,  at  these  suboptimial 
solutions,  an  “equal-ripple”  is  emerging. 

Conjecture:  SIMO  matching  solutions  are  characterized  by  an 
“equal-ripple”  gain  for  d  <  cx) 
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Gj{je>)  Gj{je>) 


ThreeWires  1 00to700MHz  matched  over  *S'( (4,4) 


Figure  16:  d  =  A  matching  over  100-700  MHz. 

ThreeWires  1 00to700MHz  matched  over  *S'( (4,4) 


Figure  17:  <7  =  4  matching  over  140-700  MHz. 
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Gj{je>)  Gj{je>) 


ThreeWires  1 00to700MHz  matched  over  *S'( /?)g  (4,3) 


Figure  18:  d  =  3  matching  over  140-700  MHz. 

ThreeWires  1 00to700MHz  matched  over  *S'( (4,2) 


Figure  19:  d  =  2  matching  over  140-700  MHz. 
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Gj{je>)  Gj{je>) 


ThreeWires  1 00to700MHz  matched  over  *S'( /?)g  (4,1) 


Figure  20:  d  =  1  matching  over  140-700  MHz. 

ThreeWires  1 00to700MHz  matched  over  *S'( (4,0) 


Figure  21:  <7=0  matching  over  140-700  MHz. 
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Table  4  tabulates  the  matching  performance  as  a  function  of  degree.  The 
rolling  off  of  the  matching  with  increasing  degree  and  the  flatness  of  the  gain 
support  the  conjecture  that  -1  dB  is  a  reasonable  matching  bound.  Indeed, 
it  is  the  flatness  of  the  gain  that  is  characteristic  of  the  theory  [21]. 
This  flat  gain  leads  to  the  excellent  research  question:  How  much  of  the  H°° 
theory  carries  over  to  the  SIMO  matching  problem? 


Table  4:  Matching  with  d)  over  140-700  MHz. 


Gain  (dB)  VSWR 


0  -1.582  3.4696 

1  -1.453  3.2849 

2  -1.368  3.1651 

3  -1.025  2.6935 

4  -1.033  2.7043 


What  is  also  intriguing  about  Table  4  is  that  the  greatest  improvement 
in  matching  is  obtained  by  the  non-reactive  4-port  for  <7  =  0.  The  associated 
scattering  matrix  is 

■  0.5748  -0.1638  -0.3528  -0.7199  ' 

0.4973  0.8192  0.2754  0.0757 

“  0.4024  -0.5057  0.7604  0.0637  ' 

0.5103  -0.2150  -0.4706  0.6870 

This  scattering  matrix  is  not  reciprocal  but  exhibits  a  non-zero  gyrator  rank 
discussed  in  Appendix  B: 

2  =  fi-iA]  =  ^rank(S'  -  S'^). 

That  is,  a  non-reactive  4-port  must  contain  at  least  2  gyrators.  Therefore, 
an  excellent  research  question  asks:  What  multiport  matching  is  possible  with 
reciprocal  (gyrator- free)  multiports? 
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6  Channelized  Multiports 

The  excellent  matching  obtained  by  the  general  multiports  (4,  d)  raises 
an  interesting  question:  are  these  multiports  masking  the  simpler  matching 
topology  of  channelized  2-ports?  Figure  22  is  a  schematic  of  the  channelized 
2-ports. 


Sjip) 


Figure  22:  Matching  with  channelized  2-ports. 

The  matching  4-port  is  obtained  as  the  compression  of  a  6-port.  Figure  23 
is  a  schematic  of  the  6-port  that  ties  Ports  1,  3,  and  5  together. 


Figure  23:  Simple  channelized  2-ports. 


The  6-port  has  admittance  matrix 


Ya®Yb®Yc  = 


VA,!!  VAA^  0 

yA,2l  yA,22  0 

0  0  yB, 11 

0  0  yB, 21 

0  0  0 
0  0  0 


0  0  0 
0  0  0 
yB, 12  0  0 

2/5,22  0  0 

0  2/C, 11  yc,i2 

0  2/C,21  yc,22 
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Using  the  labeling  in  Figure  23,  the  current  vector  for  the  4-port  is  the 
compression  of  the  current  vector  for  the  6-port: 


ii  ' 

■  1 

0 

1 

0 

1 

0  ■ 

i2 

0 

1 

0 

0 

0 

0 

is 

0 

0 

0 

1 

0 

0 

.  *4  . 

0 

0 

0 

0 

0 

1 

iA 

12 

is 

H 

ic 

ie 


Let  E  denote  the  transpose  of  this  compression  matrix  for  the  current  vectors. 
Matrix  E  is  the  expansion  matrix  for  the  voltage  vectors: 


Va  ' 

o 

o 

o 

V2 

0  10  0 

Vi 

vb 

10  0  0 

V2 

Va 

0  0  10 

Vs 

Vc 

10  0  0 

Va 

V6 

0  0  0  1 

The  admittance  matrix  Y  for  the  4-port  is  the  full  6-port  admittance  matrix 
multiplied  on  the  left  by  the  current  compression  matrix  and  on  the  right  by 
the  voltage  expansion  matrix: 


Y 


E^  {Ya®Yb®Yc}E 


yA,ii  + 1/5,11  +  yc, 11 
yA, 21 

1/5,21 

yc, 21 


yA,12  1/5,12  yc, 12 
yA,22  0  0 

0  1/5,22  0 

0  0  yc,22 
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The  preliminary  results  matching  with  this  channelized  multiport  are 
negative  but  intriguing.  Figure  24  displays  the  matching  that  shunt  inductors 
obtained  over  140-700  MHz.  The  blue  line  is  the  gain  of  the  direct  connection. 
The  green  line  is  the  gain  achieved  by  the  channelized  matching  circuit. 
Although  the  deep  null  between  140  and  300  MHz  is  hlled,  the  antenna  is 
still  poorly  matched. 


ThreeWires  100to7CX)IS/IFfe  matched,  over  channelized  ladders 


Figure  24:  Channelized  2-ports  (shunt  inductance) 

Searching  over  several  classes  of  lumped  ladders  for  better  performance 
was  fruitless.  This  empty  search  offers  two  lessons.  First,  time  and  effort 
can  be  wasted  without  good  matching  bounds.  Such  bounds  can  show  a 
design  is  simply  not  feasible  and  a  fruitless  search  is  avoided.  Second,  the 
non-reactive  matching  obtained  by  the  degree  d  =  0  multiports  shows  that 
simply  throwing  lots  of  lumped  elements  at  the  problem  is  futile — if  only 
channelized  matching  is  used.  Rather,  some  cross  coupling  must  be  present. 
Accordingly,  the  topology  was  slightly  generalized. 
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Figure  25  shows  a  power  trivider  consisting  of  channelized  2-ports  with 
admittance  matrices  Yb,  Yq,  and  Yd  that  all  terminate  in  a  common  2-port 
with  admittance  matrix  Ya- 


Figure  25:  Power  trivider. 

For  matching,  the  2-ports  are  ideal  transmission  lines.  A  uniform,  lossless 
transmission  line  of  characteristic  impedance  Zc  and  commensurate  length  i 
is  called  a  unit  element  (UE)  and  has  a  chain  matrix  [4,  Eq.  8.1] 

,  .  _  r  cosh(rp)  ZcSmh{Tp) 

^  |/cSinh(rp)  cosh(rp)  ’ 

where  r  =  c~^i  is  the  commensurate  one-way  delay  determined  by  the  speed 
of  propagation  c.  Figure  26  shows  that  this  transmission-line  topology  does 
slightly  better  than  no  matching.  The  blue  line  is  the  gain  of  the  unmatched 
gain  of  Figure  14.  The  green  line  is  the  gain  provided  by  the  matching  circuit. 
The  box  lists  the  impedances  ZcS  in  ohms  and  the  length  in  millimeters. 

Just  as  in  the  multiport  matching  of  Section  5,  better  performance  is 
obtained  by  reducing  the  frequency  range.  Figure  27  shows  pulls  in  both 
the  high  and  low  ends  to  obtain  a  3.1  VSWR  over  300-600  MHz.  The  high 
performance  obtained  at  the  600  MHz  suggests  that  the  high-frequency  range 
can  be  extended  with  a  small  increase  in  the  VSWR.  Likewise,  the  increasing 
gain  at  the  low  end  also  suggests  that  the  peak  just  below  300  MHz  can  also 
be  exploited. 


TlTree\Mres  1 00tx>700MHz:  4  ixarisrnissicri  lines 


Figure  26:  Power  trivider  matching  140-700  MHz. 


TlTree\Mres  1 00tx>700MHz:  4  ixarisrnissicri  lines 


MFfe;  II  C7^|_^= -6.533^ -1.333  dB 

Figure  27:  Power  trivider  matching  300-600  MHz. 
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Figure  28  verifies  the  observations  of  the  preceding  plots  by  extending  the 
matching  over  270-700  MHz.  The  VSWR  increases  from  3.1  to  3.5.  What 
this  example  does  show  is  that  an  explicit  matching  4-port  can  deliver  some 
performance.  Nevertheless,  this  approach  is  still  incremental  in  the  sense 
that  cross  coupling  is  not  designed  into  this  topology. 


ThreeWires  1  OOtoTOOMHz;  4  transmission  lines 


Figure  28:  Power  trivider  matching  270-700  MHz. 
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7  Exploiting  Multiport  Coupling 

This  multiport  antenna  is  an  excellent  demonstration  of  the  design  potential 
of  multiport  matching.  Exploiting  this  “lucky  hit”  has  immediate  payoffs 
for  theory,  multiport  design,  and  optimization.  Rather  than  list  a  series  of 
“research  directions”  that  constrain  cogitation,  the  following  questions  are 
offered  in  the  spirit  of  using  this  remarkable  multiport  matching  result  for 
several  points-of-departure. 

7.1  What  are  the  theoretical  limits  of  multiport  match¬ 
ing? 

This  general  question  can  be  focused  to  ask:  Given  a  multiport  load,  can 
the  bounds  or  the  Fano  bounds  be  computed  for  the  SIMO  matching 
problem?  Answering  either  theoretical  question  has  an  immediate  practical 
payoff  because  unworkable  designs  can  be  eliminated  without  wasting  time 
on  fruitless  searches  [29]. 

The  Fano  bounds  are  basic  to  single-port  matching.  A  multiport  gener¬ 
alization  was  obtained  in  1984  by  Wang  and  Chen  [40] .  Since  then,  the  Fano 
literature  shows  limited  development  for  multiports  and  is  still  restricted  to 
rational-function  models. 

The  bounds  are  basic  to  optimization  over  analytic  functions,  or 
equivalently,  to  matching  circuits.  Characteristic  of  the  bounds  is  the 
“flatness”  of  the  objective  function  at  a  local  optima;  Helton  observed 
that  the  flatness  and  the  winding  number  characterize  local  optima 
[21,  Theorem  9.3.1].  In  this  context,  the  observed  flatness  in  Figure  17 
provides  graphical  evidence  of  a  theoretically  optimal  match. 

Further  consideration  of  the  SIMO  question  also  reveals  new  approaches 
to  constrained  optimization.  Section  5  sets  out  the  SIMO  matching 
problem  for  A^-port  loads  Sl  d.s 

sup{||Gr(sG, -S',  S'ijlln-oo  :  S'  e  6}, 

where  ©  is  the  class  of  designable  (A^  -|-  l)-ports.  If  the  circuit  designer  is 
searching  for  optimal  matching  from  lumped,  lossless  multiports  of  degree 
not  exceeding  d,  the  set  inclusion 

6  C  U+{N  +  l,d), 
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implies  the  inequality: 


snp{\\GT{sG,S,SL)\\n,-oo  :  S  e  &} 

<  sup{||Gr(-SG,  S,  S'l)||j7-oo  :  S  e  U^{N  +  1,  d)}. 

Because  f/’''(iV  + 1,  d)  is  a  subset  of  several  techniques  can  com¬ 
pute  upper  bounds  on  this  gain.  One  upper  bound  follows  from  maximizing 
the  gain  for  over  the  unit  ball  in  H°°: 


SIMO  Matching:  Given  the  generator’s  reflectance 
sg  and  an  iV-port  load  Sl{p),  maximize  the  transducer 
power  gain 

maximize:  ||Gr(sG, -S',  S'l)||q  _oo, 
where  S  is  an  function  of  norm  ||5'||oo  <  1- 


Helton  solved  this  multiport  matching  problem  in  his  1981  paper  [20, 
Test  IV].  Helton’s  multiobjective  optimization  program  can  also  solve  this 
problem  [22].  The  idea  is  to  compute  the  Pareto  front  constructed  of  the 
gain  and  multiport  norm.  This  Pareto  front  subsumes  the  preceding 
SIMO  matching  as  a  special  case.  Consequently,  comparing  both  methods 
applied  to  multiport  matching  is  an  excellent  research  topic. 

More  attuned  to  Lagrange  multipliers  is  the  following  constrained  opti¬ 
mization  problem: 


maximize:  ||Gt(sg, -S',  S'ijUn -oo 

subject  to  the  constraints  that  S  is  an  function  of  unit  norm: 

||/-^^^||oo  =  l. 

A  reasonable  conjecture  is  that  the  gradient  alignment  conditions  of  the 
Lagrange  multipliers  have  an  analog  in  optimization.  Regardless  of  the 
methods,  computing  best  possible  bounds  for  multiport  matching  provides  a 
fundamental  benchmark  to  assess  practical  design.  However,  there  remains 
the  hard  problems  of  actually  synthesizing  optimal  and  near-optimal  match¬ 
ing  multiports. 
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7.2  What  are  the  optimal  multiport  matching  circuits? 

Figures  17-21  show  near-optimal  matching  performance  of  the  lumped,  loss¬ 
less  4-ports.  Consequently,  the  extraction  of  a  circuit  from  these  near-optimal 
multiports  is  a  practical  problem  of  synthesis.  One  approach  is  the  partial- 
fraction  or  Foster  expansion.  Although  the  Foster  expansion  does  produce 
a  multiport  circuit,  these  circuits  are  not  practical  because  they  are  infested 
with  transforms  and  gyrators — the  value  of  this  approach  is  to  demonstrate 
feasibility  [35].  Therefore,  other  approaches  to  extract  a  practical  circuit 
from  these  lumped,  lossless  4-ports  are  credible  lines  of  research.  The  payoff 
from  these  considerations  have  immediate  applications:  antenna  matching 
[34],  ampliher  design,  and  circulator  optimization  [45]. 

Specihc  to  the  three-wire  antenna,  the  multiports  need  only  a  few  lumped 
elements  to  get  a  good  match.  In  particular,  the  multiport  of  Figure  21  has 
no  lumped  elements.  Its  associated  scattering  matrix  is  constant  and  listed 
in  Section  5.  Decompositions  of  this  constant  scattering  matrix  may  reveal 
the  basic  connections  between  the  ports  [28,  Section  7-2]  or  relate  to  the 
decoupling  algorithms  of  Volmer  [37]. 


7.3  How  effective  is  decoupling? 

Figure  1  illustrates  the  decoupling  algorithm  [37].  The  idea  is  excellent — use 
the  eigenvectors  of  the  load’s  scattering  matrix  to  decouple  the  multiport 
load  into  a  collection  of  single-port  loads.  The  key  word  in  the  question  is 
“effective”  because  problems  exist  with  synthesizing  the  decoupling  circuit 
over  wide  frequency  band.  The  basic  idea  starts  with  the  observation  that  if 
Sl  is  an  A^-port  that  diagonalizes  as 

A  =  diag(Ai,A2,...,A^)  =  77^5^17, 

where  U  is  an  N  x  N  unitary  matrix,  the  2A^-multiport 


Su  = 


0 

U  0 


is  both  lossless  and  decouples  the  load  in  cascade: 

A  =  d{Su,SL)  =  U^SLU. 

The  technical  point  of  this  technique  is  to  synthesize  the  decoupler.  Volmer’s 
approach  uses  the  2-port  directional  couplers  to  synthesize  Su  similar  to  the 
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way  Givens  rotations  can  diagonalize  a  matrix.  This  deconpling  techniqne 
is  fonnd  in  other  mnltiport  applications.  For  example,  Czawka  [10]  offers 
broadband  microstrip  nnconplers  for  mnltiport  complex  loads.  Closely  allied 
with  this  SIMO  design  is  a  six-channel  wideband  power  divider  by  Czawka 
[11].  Jarmasz  and  Martens  [24]  apply  a  variant  of  this  deconpling  approach 
to  wave  digital  hlters.  Finally,  Cameron  and  Yn  [6]  offer  several  mnltiplexer 
topologies  that  can  be  applied  to  deconpling  design.  For  example,  a  wideband 
deconpler  may  be  obtained  by  channelizing  several  narrowband  deconplers 
across  the  freqnency  band. 

7.4  What  are  the  multiport  ladder  circuits? 

The  lossless  ladders  are  basic  to  2-port  matching.  Each  ladder’s  stage  com¬ 
monly  consists  of  a  series  or  shnnt  element.  The  stages  are  cascaded  to 
make  a  ladder  of  specihed  degree.  By  increasing  the  degree  of  the  ladder, 
the  matching  performance  of  the  ladder’s  topology  (low-pass,  high-pass)  as 
a  fnnction  of  degree  can  be  benchmarked  against  the  bonnds  [1] .  There¬ 
fore,  generalization  of  these  2-port  ladders  to  mnltiport  ladders  is  a  credible 
research  topic  with  immediate  and  practical  payoffs.  As  illnstrated  in  Fig- 
nre  29,  one  mnltiport  generalization  nses  a  stage  consisting  of  shnnts  between 
ports  [31],  [8].  Jnst  as  Fnjisawa  obtained  deep  resnlts  on  characterizing  the 
mid-series  and  mid-shnnt  ladders  [4],  and  as  Fialkow  pnshed  this  program  to 
parallel  ladders  [16],  an  analogons  mnltiport  ladder  program  that  bnilds  on 
these  classic  resnlts  awaits  the  determined  researcher. 


Fignre  29:  Shnnts  between  ports. 
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7.5  How  critical  is  multiport  matching  for  MIMO? 

Figure  30  illustrates  multiport  matching  in  Multiple-Input  Multiple-Output 
(MIMO)  systems.  Although  maximum  power  transfer  is  the  standard  match¬ 
ing  objective,  MIMO  system  performance  may  be  enhanced  more  by  decou¬ 
pling  the  antenna  rather  than  power  matching. 


Matching  Channel  Matching 


Figure  30:  MIMO  matching  networks  [25]  ©  IEEE  2006. 

Fei,  Fan,  and  Thompson  offer  a  “multi-handed”  discussion  of  MIMO 
matching  [17]: 

. . .  integration  of  MIMO  technique  into  compact  devices  is  re¬ 
stricted  if  the  antenna  spacing  is  below  half  a  wavelength.  This 
because  strong  mutual  coupling  (MC)  between  closely  spaced  an¬ 
tenna  elements  results  in  changes  in  antenna  patterns  (antenna 
correlation)  and  loss  of  antenna  efficiency  . . .  MC  is  claimed  as 
a  detriment  to  MIMO  systems.  However,  MC  can  also  be  a 
positive  factor  to  increase  the  MIMO  performance  under  some 
circumstances. 

These  contrapuntal  arguments  were  addressed  by  Wallace  and  Jensen 
[39]  showing  MIMO  capacity  matching  is  optimized  by  conjugate  matching. 
However,  Fei,  Fan,  and  Thompson  also  observe  that  [17]: 

...  it  is  not  feasible  to  integrate  this  solution  [conjugate-match] 
into  MIMO  systems  as  multi-loads  has  to  be  introduced  to  each 
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receive  port . . .  Although  the  design  of  multiport-conjugate  match¬ 
ing  network  has  been  reported  in  . . .  the  authors  do  not  know  of 
any  experimental  results  being  presented  in  the  literature.  As  a 
result,  the  single-port  match  is  still  an  attractive  if  suboptimal 
solution  .... 

The  upshot  is  that  single-port  matching  is  a  credible  benchmark  for  proposed 
multiport  matching — any  MIMO  matching  scheme  must  beat  this  subopti¬ 
mal  performance.  The  upper  bound  for  MIMO  matching  is  the  conjugate 
matching  [39].  Consequently,  there  are  excellent  multiport  design  opportu¬ 
nities  in  wideband  MIMO.  Finally,  other  multiport  design  issues  arise  when 
transmit  and  receive  technologies  are  integrated  into  small  systems.  Quoting 
from  Morris  and  Jensen  [26]: 

As  adoption  of  MIMO  technology  increases,  there  will  be  in¬ 
creased  desire  to  integrate  multiple  receiver  front  ends  on  a  single 
chip,  particularly  for  mobile  equipment.  As  this  integration  oc¬ 
curs,  circuit  level  signal  coupling  will  increase,  potentially  leading 
to  altered  signal  correlation  characteristics  and  signal-to-noise  ra¬ 
tio  (SNR)  at  the  front  end  amplifier  outputs. 

Figure  31  shows  a  multiport  embedded  in  an  active  MIMO  system.  Mor¬ 
ris  and  Jensen  show  that  multiport  matching  can  reduce  the  electromagnetic 
signal  coupling  in  the  radio  receiver  front  end.  Because  the  amplifiers  will  in¬ 
troduce  gain,  noise,  and  stability  problems,  scalar-valued  multiport  optimiza¬ 
tion  turns  into  vector-valued  multiport  optimization  and  the  computation  of 
the  associated  Pareto  fronts. 


Figure  31:  MIMO  receiver  and  matching  network  [26]  ©  IEEE  2005. 
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7.6  Multiport  Antenna  System  Optimization 

Up  to  this  last  example,  multiport  matching  has  been  focused  on  matching 
the  3-port  load  as  a  black  box  rather  than  as  an  antenna.  This  subsection 
considers  the  performance  of  the  antenna  system  both  “in  the  wires”  and 
“over-the-air.”  That  is,  the  matching  multiport  and  the  antenna’s  geometry 
are  simultaneously  twiddled  to  optimize  system  performance.  For  example, 
Czawka  and  Garbaruk  [12]  offer  a  design  framework  to  control  the  radiation 
pattern  from  a  multiport  antenna  by  using  a  lossless  multiport  to  control 
phase  of  the  feeds.  A  fast  and  generalized  antenna  code  is  critical  to  the 
optimization  of  such  an  antenna  system  [7]. 

Figure  32  illustrates  such  a  system  optimization  and  blurs  the  distinction 
between  the  matching  circuit  and  the  antenna.  The  antenna  is  still  the 
three-wire  antenna.  The  antenna  system  is  driven  by  wire  fed  through 
its  matching  2-port  Si.  The  short  wire  7^2  and  the  long  wire  7^3  are  fed 
by  electromagnetic  coupling.  These  two  wires  are  loaded  at  their  feed  points 
with  2-ports  with  scattering  matrices  S2  and  S3.  The  design  variables  are  the 
three  2-ports — and  the  parameters  of  the  wires:  their  lengths  and  spacings. 
The  design  problem  is  to  simultaneously  match  while  sweeping  over  this 
family  of  three-wire  antennas. 


Figure  32:  Antenna  system  (single  feed,  port  loading). 

One  optimization  problem  is  solved  by  computing  a  Pareto  front  that 
shows  the  tradeoff  between  gain  and  bandwidth.  A  more  ambitious  opti- 
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mization  problem  is  solved  by  computing  a  Pareto  surface  that  shows  the 
tradeoffs  between  gain,  bandwidth,  and  volume.  If  a  radiation  pattern  is 
also  specified,  the  optimization  problem  is  solved  by  computing  the  Pareto 
surface  that  shows  the  tradeoffs  between  gain,  bandwidth,  and  volume — 
under  the  constraint  that  the  resulting  radiation  pattern  is  within  tolerance. 
A  handy  “rule-of-thumb”  in  optimization  theory  is  that  any  constraint  is 
really  an  objective  fnnction.  Conseqnently,  the  hnal  antenna  design  prob¬ 
lem  is  computation  of  the  Pareto  snrface  showing  the  tradeoffs  between  gain, 
bandwidth,  volnme,  and  qnality  of  pattern  match.  Clearly,  a  fast  mnlti- 
physics  code  is  needed  to  handle  the  Inmped  elements  in  the  circnit  and  the 
distribnted  electromagnetic  performance  of  the  antenna. 
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A  Characteristic  Modes  of  the  Three-Wire 
Antenna 

Because  the  coupling  between  the  ports  of  the  three-wire  antenna  gives  its 
wideband  performance,  the  discussion  of  Section  3  is  augmented  by  this  ap¬ 
pendix.  The  characteristic  modes  reveal  the  resonances  and  couplings  of  the 
three  wires.  Figure  A-1  re-plots  the  three  wires  for  the  reader’s  convenience. 


Figure  A-1:  Center-fed  Dipoles:  1/4,  1/2,  and  1  meter  in  length 

Figure  A-2  displays  the  impedance  matrix  determined  by  the  segmenting 
of  the  wires  as  listed  in  Table  3.  The  matrix  breaks  into  three  blocks  cor¬ 
responding  to  the  three  wires.  The  blocks  on  the  main  diagonal  show  the 
coupling  between  the  segments  on  a  wire.  The  off-diagonal  blocks  show  the 
impedance  coupling  between  the  wires.  Write  the  impedance  matrix  using 
its  real  and  imaginary  parts: 

Z  =  R  +  ]X. 

The  resistance  and  reactance  are  real  matrices,  with  i?  >  0.  Make  the  as¬ 
sumption  that  the  segmentation  is  sufficiently  homogenous  so  that  these  real 
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Impedance  Matrix:  ThreeWires  150MHz 
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Figure  A-2:  Impedance  matrix  at  150  MHz. 

matrices  are  also  symmetric.  The  characteristic  modes  are  the  eigenvalues  of 
[19]: 

XJ  =  RJA]  A  =  diag(Ai, . . . ,  Aat), 

where  J  is  matrix  of  eigencurrents.  Because  X  and  R  are  real  and  symmetric, 
A  must  be  real.  If  R  is  invertible  {R  >  0),  replace  J  =  R~^/'^W  to  get 

R^/^XR-^/^W  =  WA. 
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This  ordinary  eigenvalue  problem  admits  a  solution  W  that  is  real  and  or¬ 
thonormal:  /  =  W'^W.  Consequently,  the  eigencurrents  J  are  orthonormal 
with  respect  to  the  resistance  matrix: 

J^RJ  =  I. 

A  physical  meaning  of  the  characteristic  modes  follows  from  the  impedance 
matrix 

ZJ  =  RJ{I  +jA). 

Resonance  is  said  to  occur  when  a  characteristic  mode  vanishes  [23]:  =  0. 

If  \n  >  0,  its  eigencurrent  jn  is  said  to  make  an  inductive  contribution.  If  A„  < 
0,  its  eigencurrent  j„  is  said  to  make  a  capacitive  contribution.  Properties  of 
the  eigencurrents  are  offered  by  Austin  and  Murray  [3] : 

An  important  property  of  these  current  modes  is  that  they  are 
orthogonal  over  the  surface,  as  are  the  helds  that  they  produce 
over  the  sphere  at  inhnity,  and  so  they  radiate  power  indepen¬ 
dently  of  one  another.  The  magnitude  of  A„  denotes  the  relative 
signihcance  of  each  mode,  in  terms  of  the  ratio  of  stored  energy  to 
radiated  power.  The  closer  |A„|  is  to  zero,  and  hence  to  resonance, 
the  more  substantial  is  its  contribution  to  the  overall  radiation. 

If  these  modes  can  be  individually  excited  by  an  adroit  placement  of  feed 
points,  design  of  antennas  and  shaping  of  radiated  helds  follows  from  the 
orthogonality  of  the  eigencurrents.  The  following  plots  illustrate  these  con¬ 
cepts. 

Figure  A-3  plots  the  two  most  radiating  eigencurrents  at  150  MHz.  Be¬ 
cause  the  wire  segments  are  the  same  lengths  and  the  indices  run  along  the 
wires,  the  plot  shows  the  current  along  the  wires.  The  half-wave  resonance 
of  the  long  wire  with  a  characteristic  mode  of  A  =  0.61187  (blue  curve). 
Most  of  the  current  resides  on  the  long  wire  but  there  is  current  on  both  the 
shorter  wires.  Comparison  with  Figure  6  reveals  150  MHz  is  sightly  above 
resonance  of  the  long  wire.  Consequently,  the  characteristic  mode  is  slightly 
positive  rather  than  zero.  If  a  feed  point  is  placed  at  the  center  of  the  long 
wire,  this  mode  could  be  excited  at  150  MHz  and  radiate  most  of  its  power. 
Figure  A-3  also  plots  the  next  most  radiative  eigencurrent  (A  =  —143)  that 
corresponds  to  the  full-wave  resonance  of  the  long  wire.  If  feed  points  are 
placed  at  the  extrema  of  this  eigencurrent,  this  mode  will  be  excited  (with 
opposite  phasing)  at  150  MHz.  However,  the  radiation  from  this  eigencurrent 
with  A  =  —143  is  relatively  small. 
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Figure  A-3:  Eigencurrents  at  150  MHz. 
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Figure  A-4  plots  two  significant  eigencurrents  at  300  MHz.  The  medium- 
length  wire  does  show  half-wave  resonance  (A  =  0.65159)  and  shows  coupling 
into  the  other  wires.  However,  the  most  efficient  resonance  (A  =  0.47758) 
is  the  full-wave  resonance  of  the  long  wire  (blue  line).  Thus,  300  MHz  is 
slightly  above  resonance  for  both  the  long  and  medium-length  wires.  A  feed 
point  placed  at  the  center  of  the  medium-length  wire  should  excite  this  mode. 
Two  feed  points  on  the  long  wire  and  oppositely  phased  should  drive  the  full- 
wave  resonance.  Weighting  and  phasing  these  two  eigencurrents  allows  some 
beamforming  in  the  far  field. 


ThreeWires  300MHz.ZMM:  Eigencurrents  at  300  MHz 
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Figure  A-4:  Eigencurrents  at  300  MHz. 
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Figure  A-5  plots  three  significant  eigencurrents  to  600  MHz:  the  half 
wave  on  the  short  wire,  the  full  wave  on  medium-length  wire,  and  the  two 
waves  on  the  long  wire.  This  plot  reveals  some  interesting  feed-point  design 
problems  with  respect  to  exciting  individual  modes.  A  single  feed  centered 
on  the  short  wire  should  only  excite  its  half  wave  (A  =  0.45367).  Two  feeds 
on  the  medium-length  wire  (A  =  0.44225)  produce  the  full-wave  on  that 
wire  and  two  waves  on  the  long  wire.  Likewise,  four  feeds  on  the  long  wire 
should  excite  the  two-wave  current  (A  =  0.21943).  Thus,  a  total  of  seven 
feeds  should  allow  control  of  these  three  eigencurrents  and  the  associated 
orthogonal  decomposition  of  the  far  field. 


ThreeWires  600MHz.ZMM:  Eigencurrents  at  600  MHz 
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Figure  A-5:  Eigencurrents  at  600  MHz. 
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Figure  A-6  shows  significant  eigencurrents  at  450  MHz.  The  plot  reveals 
that  the  three-wire  antenna  is  radiating  from  the  long  wire  (A  =  0.36559) 
with  3/2  wave.  The  extrema  of  this  wave  should  set  the  location  of  the  feeds 
on  the  long  wire.  However,  the  A  values  of  the  other  eigencurrents  show 
radiation  at  450  MHz  could  be  accomplished  from  several  locations. 


ThreeWires  450MHz.ZMM:  Eigencurrents  at  450  MHz 
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Figure  A-6:  Eigencurrents  at  450  MHz. 

In  summary,  these  plots  demonstrate  that  the  three-wire  antenna  has  sev¬ 
eral  radiation  modes,  shows  the  coupling  between  the  wires,  locates  credible 
feed  points  and  phasing  for  this  antenna  system.  Thus,  the  design  insight 
provided  by  the  characteristic  modes  by  locating  these  feed  points  and  phas¬ 
ing  is  orders-of-magnitude  more  efficient  that  the  “ant-like”  combinatorial 
algorithms  that  search  all  segments  and  phases  for  such  control  points.  Fi¬ 
nally,  the  design  insight  provided  by  the  characteristic  modes  also  point  to 
the  2-port  loading  led  to  Figure  32  and  points  to  end-to-end  loading  with 
low-pass  hlters  for  better  performance  at  the  lower  frequencies. 
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B  Lumped,  Lossless  A^-Ports 


Let  {N,  d)  denote  the  class  of  lumped,  lossless  iV-ports  of  degree  d.  Phys¬ 
ically,  U~^  [N,  d)  represents  all  the  iV-ports  that  can  be  synthesized  with  at 
most  d  lumped  elements — inductors  and  capacitors — connected  only  by  wires, 
transformers,  and  gyrators.  Excellent  wideband  matching  to  the  three-wire 
antenna  was  obtained  from  the  4-ports  with  degree  d  >  3.  This  appendix 
reviews  U~^  (TV,  d)  to  analyze  these  matching  results. 

Mathematically,  U~^{N,  d)  consists  of  all  iV  x  iV  scattering  matrices  S{p) 
that  are  real,  rational,  analytic  and  bounded  on  the  open  right  half  of  the 
complex  plane  C+,  unitary  on  the  imaginary  axis 

S{juj)^  S{juj)  =  In, 
and  have  Smith-McMillan  degree 


degsM[>^(p)]  <  d. 


The  initial  results  on  the  Belevitch’s  Theorem  characterizes  the  scattering 
matrices  for  the  lumped,  lossless  2-ports  of  U~^{2,  d). 

A  lumped,  passive  2-port  contains  the  lumped  elements — inductors,  ca¬ 
pacitors,  resistors — and  the  means  of  connections — wires,  transformers,  and 
gyrators.  The  corresponding  scattering  matrix  [28],  [42] 


S{p) 


Snip)  Si2{p) 

S2l{p)  S22{p) 


is  a  rational  function  that  is  analytic  on  the  open  right  half  plane  C+,  satishes 
the  real  condition 

Sip)=W)  (peC+), 

and  is  a  contraction 


S{p)^S{p)  < 


1 

0 


A  lumped,  lossless  2-port  omits  the  resistors.  As  a  consequence,  its  scat¬ 
tering  matrix  S{p)  is  a  real  rational  function  that  is  analytic  on  C+,  a  con¬ 
traction  on  C+,  and  unitary  on  the  boundary  of  C+: 

S{juj)^S{juj)  = 


1  0 
0  1 
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A  lumped,  passive,  reeiprocal  2-port  omits  the  gyrators.  The  associated 
scattering  matrix  S{p)  is  a  real  rational  function  that  is  analytic  on  C+,  and 
a  contraction  on  C+  that  is  also  symmetric: 

S(p)  =  S(pf  (p€C+). 

Belevitch’s  Theorem  asserts  the  existence  of  a  scattering  matrix  S(p)  for  each 
2-port  and  obtains  the  converse  of  these  statements. 


Belevitch’s  Theorem  [44],  [4,  pages  83-86]  A  lumped,  lossless 
2-port  admits  a  scattering  matrix 

h(p)  fip)  1 

e/*(p)  -^K{p)  \  ’ 

where  h*(p)  =  h{—p)  and  (/,  g,  h)  is  a  Belevitch  triple: 

B-1  /(p),  g{p),  and  h{p)  are  real  polynomials;  e  =  ±1, 

B-2  g{p)  is  strict  Hurwitz  (no  zeros  in  the  closed  right  half  plane). 
B-3  g^{p)g{p)  =  Mp)f{p)  +  K{p)h{p)  for  all  p  e  C. 


Conversely,  any  such  S{p)  is  the  scattering  matrix  for  some  lumped, 
lossless  2-port.  Moreover,  if  the  2-port  is  reciprocal 


S{p)  = 


1 

9{p) 


Kp)  fip) 

f{p)  -eK{p) 


•> 


where  either  /(p)  is  even  and  e  =  —1,  or  /(p)  is  odd  and  e  = 
1.  Conversely,  any  such  S{p)  is  the  scattering  matrix  for  some 
lumped,  lossless,  reciprocal  2-port. 


The  number  of  reactive  lumped  elements — the  inductors  and  capacitors — 
is  linked  to  the  degree  of  the  2-port. 

Definition  1  page  91]:  Let  the  K  distinet  poles  of  S{p)  he  denoted  as 
Pk-  Let  deg(S'(p);pfc)  denote  the  largest  order  to  whieh  p  =  Pk  appears  in  any 
minor  of  S{p).  The  Smith-McMillan  degree  of  S{p)  is 

K 

degsM[*^(p)]  :=  5^deg(A(p);pfc). 

k=l 
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Belevitch’s  Theorem  characterizes  U~^{2,d).  The  scattering  matrices  in 
U~^{2,d)  constitute  the  “state-space”  curve  in  Figure  12.  Although  Bele¬ 
vitch’s  Theorem  admits  some  direct  extensions  [27],  the  state-space  repre¬ 
sentation  is  more  useful  for  A^-ports.  Figure  B-1  illustrates  the  state-space 
representation  of  a  passive,  lumped  A^-port. 


Figure  B-1:  State-space  representation  of  a  passive,  lumped  A^-port. 

The  hgure  shows  that  by  pulling  the  d  reactive  elements  and  the  r  lossy 
elements  into  the  augmented  load  Sl{p),  what  is  left  is  an  multiport  consisting 
of  the  wires,  transformers,  and  gyrators  with  scattering  matrix  Sa-  Because 
Sa  represents  a  multiport  containing  only  wires,  transformers,  and  gyrators, 
Sa  is  a  constant  matrix.  Moreover,  because  the  multiport  containing  only 
wires,  transformers,  and  gyrators  is  lossless,  Sa  is  a  constant  unitary  matrix. 

The  A^-port  is  obtained  by  looking  into  Ports  1,  2,  . . . ,  A^  of  the  aug¬ 
mented  scattering  matrix  Sa  while  its  remaining  ports  are  terminated  in  the 
augmented  load  Sl{p)-  That  is,  S{p)  is  the  image  of  the  augmented  load 
viewed  through  the  augmented  scattering  matrix  Sa-  Theorem  1  gives  the 
precise  statement  of  this  state-space  representation. 
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Theorem  1  (State-Space  [42,  pages  90—93])  Every  lumped,  passive,  ca¬ 
sual,  time-invariant  N -port  admits  a  scattering  matrix  S{p)  and  conversely. 
If  S{p)  has  Smith- McMillan  degree  d  and  normal  rankr,  defined  as  [42,  page 
91]: 

r[S{p)]  :=  rank[/Ar  -  S'(-p)^S'(p)], 

S{p)  admits  the  state-space  representation; 


S{p)  :=  Sa,ll  +  Sa,l2SL{p){Id+r  "  ^a,22^L (p) ) , 


where  the  augmented  load  is 


Slip) 


q  X  0  0 

0  —q  X  I i\f^  0 

0  0  Ox/,. 


and  Nl  -\-  Nc  =  d.  The  augmented  scattering  matrix  is 


Sa 


Sa,ll 
Sa, 21 
N 


Sa, 12 
Sa,22 
d  +  r 


N 

d  +  r 


which  is  a  constant,  real,  orthogonal  matrix.  If  the  N -port  is  lossless,  then 
r  =  0.  If  the  N -port  is  reciprocal,  then  Sa  is  symmetric,  Sfi  =  Sa. 


The  power  of  the  state-space  representation  is  that  the  scattering  matrix 
encodes  the  structure  of  the  iV-port  [42,  page  91],  [28]: 

•  Any  A^-port  that  has  S{p)  as  its  scattering  matrix  contains  at  least 
r[5'(p)]  resistors.  Moreover,  of  all  the  iV-ports  that  have  S{p)  as  their 
scattering  matrix,  at  least  one  has  exactly  r[5'(p)]  resistors. 

•  Let  Nl  and  Nc  denote  the  number  of  inductors  and  capacitors  in  an 
7V-port.  For  every  iV-port  with  scattering  matrix  S{p),  Nl  +  Nc  > 
deggM[>S'(p)].  Moreover,  of  all  the  iV-ports  that  have  S{p)  as  their  scat¬ 
tering  matrix,  at  least  one  iV-port  has  exactly  this  many  inductors  and 
capacitors. 

•  The  number  of  gyrators  in  any  iV-port  must  always  exceed  the  gyrator 
rank: 

9[Sip)]  ■■=  ^rank[A(p)  -  ^(p)^]. 
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Oono  [30]  established  that  exactly  5f[S'(p)]  gyrators  are  needed  to  syn¬ 
thesize  the  iV-port.  Consequently,  when  5f[S'(p)]  =  0  or  when  S{p)  is 
symmetric,  the  iV-port  is  reciprocal  and  contains  no  gyrators.  Phase  2 
of  this  series  will  report  that  excellent  wideband  matching  is  still  possi¬ 
ble  using  reciprocal  multiports.  Specihcally,  good  matching  is  obtained 
using  only  capacitors — no  inductors  and  no  gyrators. 

•  Anderson  [2]  obtains  the  general  result  regarding  iV-port  synthesis:  an 
iV-port  synthesis  is  possible  with  deggjy[[S'(p)]  reactive  elements  and 
5'[*S'(p)]  gyrators. 

Section  5  applies  this  state-space  representation  for  wideband  matching  of 
the  three-wire  antenna  by  specializing  to  the  d). 
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significant  research  topics.  Another  direction  aims  at  the  matching  networks  because  of  the  gaps  in  theory  and  practical  design. 
Specific  and  detailed  formulation  of  these  research  opportunities  is  made  explicit  in  this  report. 
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multiport  matching  wideband  frequency  channelized  multiports 

port  decoupling  three-wire  antenna  Pareto  surface 


|16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

J.  C.  Allen 

19B.  TELEPHONE  NUMBER  (Include  area  code) 

u 

u 

u 

uu 

64 

1619')  553-6566 

Standard  Form  298  (Rev.  8/98) 


Prescribed  by  ANSI  Std.  Z39.18 


Approved  for  public  release;  distribution  is  unlimited. 


